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PREFACE 


The work described in this report was performed by the Applied 
Mechanics Division of the Jet Propulsion Laboratory under the cognizance 
of the Mariner Jupiter /Saturn 1977 Project. The work was originally pre- 
pared as a paper which was presented at the Winter Meeting of the American 
Nuclear Society held in San Francisco, California, on November 16-21, 
1975. 
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ABSTRACT 


The computing software that was used to perform the charged particle 
radiation transport analysis and shielding design for the Mariner Jupiter/ 
Saturn IS?”^ spacecraft is described. Electron fluences, energy spectra and 
dose rates obtained with this software are presented and conipared with 
independent computer calculations. 
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EEPKODUCIBiLriT OF 'JTIE 
ORIGINAL FAG r 13 pnnn 

I. INTRODUCTION 

In situ measurements of the Jovian trapped radiation were made by 
Pioneer 10 in December 1973 and by Pioneer il in December 1974. These 
measurements revealed a potentially hazardous environment for MJS'77 
spacecraft electronics and some surfaces. As a result, the Mariner 
Jupiter /Saturn 1977 project carried out an intensive study of the charged 
particle environment derived from those measurements and the effects of 
that environrr:ent on both spacecraft and mission design (Ref. 1), This 
report covers one aspect of that study: the radiation transport analysis 
and shielding design. 

Radiation analyses were performed for virtually all of the engineering 
and science subsystems of the MJS'77 spacecraft (Fig. 1) with the computer 
programs described in this r eport. System level and subsystem level shielding 
for the MJS'77 spacecraft were based on these radiation shielding calculations. 

II. SUMMARY 

A. DATA FLOW 

A flow diagram of the main computer programs is shown in Fig. 2. 
Unshielded electron and proton radiation environments are Input to the 
SHIELD program in the form of flux and fluence energy spectra for a particu- 
lar trajectory derived from Jupiter electron and proton and solar flare pro- 
ton radiation models (Ref. 1). The SHIELD program outputs include one- 
dimensional attenuation kernels .such as dose or fluence vs absorber 
thickness for electrons, protons and secondary photons (bremsstrahlung). ^ 

This information, along with the spacecraft and/or .'•ubsystem geometry 
and mass distribution descriptions, is input to the SIGMA program (Ref. 2), 
which uses three-dimensional ray tracing to obtain dose and fluence esti- 
mates at internal spacecraft locations. The SOCODE program then uses the 
SIGMA shield sensitivity data in combination with a predetermined set of 
radiation level criteria to obtain an optimum shield configuration. 

^Dose is generally expressed in units of rads, where a rad represents 
100 ergs of energy deposited per gram of material. In this report the 
material is silicon. 
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B. 


SOFTWARE 


Several charged particle shielding codes were developed and/or 
revised during the course of thir- work. The SHIELD program, written at 
the Jet Propulsion Laboratory, subdivides materials into many small 
layers; within each layer, incident particl -s are transported using the con- 
tinuous slowing down approximation for a set of solid angle bins. The angu- 
lar distribution is revised at the surface of each layer using angular 
straggling distributions. The revised distribution is then used as the inci- 
dent data for the next layer. SHIELD uses the BETA-II program (Ref. 3) 
cross section and source spectra processors. Output includes angular 
fluxes and flux responses with slab and center-of-sphere kernels, both 
printed and punched. Results for both differential and integral fluence are 
in excellent agreement with calculations performed using the BETA-II and 
SANDYL (Ref. 4) Monte Carlo computer programs. 

The original SIGMA complex geometry program interpolated center-of- 
sphere kernels using slant path mass thicknesses for each ray used in the 
solid angle integration. This kernel agrees with Monte Carlo calculations 
for spherical geometry but underestimates radiation levels for slab geom- 
etry. A revised kernel was implemented which agrees with Monte Carlo 
calculations for both slab and spherical geometries. T.’ie revised kernel 
uses perpendicular thickness (estimated from slant paths and normal 

derivatives at material boundaries) in conjunction with a depth-dependent 
power law for angular transmission. 

The geometry description and ray tracing portions of SIGMA received 
extensive modifications, including simple input for multiple bay spacecraft, 
recognition of simple geometric shapes with multiple bounding surfaces, 
a.id implementation of CRT plotting capabilities. In addition, the require- 
ment tor explicit description of void volumes was removed.^ 

The shield s'.nsitivity option of SIGMA was exte.nded to include auto- 
matic recognition and retention of shield crossing combinations. This sensi- 
tivity information is used for designing spot shields for critical components 
and/or increasing vehicle surface thicknesses in an optimum manner. 
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III. ONE -DIMENSIONAL CALCULATIONS 


A. SLANT FATH AND CENTER-OF-SPHERE KERNELS 

A. typical one-dimensional attenuation kernel is the radiation level, 
^sphere^^^' at the center of a sphere of radius x, due to an externally inci- 
dent, cosine-distribu' ed radiation source (isotropic flux environment). This 
is the kernel used by the original SIGMA program. The kernel can be gen- 
erated by several methods and is supplied to the three-dimensional program 
as a tabulation of radiation level vs mass thickness. 


The CHARGE program (Ref, 5) is a one-dimensional code which uses 
the basic range-energy relation, modified by applying electron transmission 
factors derived from curve fit Monte Carlo data (Mar formula. Ref. 6). 

The BETA-II program uses Monte Carlo methods to predict electron 
transmission. Center-of- sphere results were approximated in slab geometry 


where x is the mass thickness of a slab, 0 is the off-normal scattering angle, 
and AD (x, cos 0) is the radiation level contribution from the forward directed 
portion (0. 9 < cos 6 < 1. 0) of the electron transmission; i. e. , only those 
particles transmitted within AH steradians of the slab normal are counted. 

Figure 3 is a comparison of a one-dimensional kernel as generated for 
the Jupiter electron environment (Ref. 1) using the CHARGE and BETA-II 
programs. Both problems assume incident isotropic flux and cosine source 
angular dependence for an aluminum spherical shell absorber with the dose 
point at the center. 

Excellent agreement between CHARGE and BETA-II was obtained 
2 2 

through 4 g/cm . Divergence beyond 4 g/cm is attribute! to the treatment, 
by CHARGE, of the high-energy ( >10 MeV) content of the Jupiter electron 
environment. First, the Mar transmission formula was based on earth 
trapped radiation energies; i. e, , <10 Mev. Second, CHARGE assumes that 
electrons slow down while travelirig :n straight lines (i. e. , straight ahead 
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approximation). This assumption leads to overly conservative flu-’iice 
estimates. 

B. SHIELD PROGRAM 

The SHIELD program was written to retain the cost effectiveness of 
CHARGE kernel generation while removing the deficien.ies for energie- 
greater than 10 MeV, SHIELD uses the BETA-II electron /photon cross 
section processing; however, kernels are generated by numerically integrat- 
ing the one-dimensional transport equation. 

Numerical methods used in the SHIELD program include (see 
Appendix A): 

(1) Subdividing material layers into many differential sublayers. 

(2) Using condensed history angular straggling distributions 
(Goudsmit - Saunderson method) for each differential sublayer. 

(3) Regrouping of electrons after each differential sublayer into a 
-.’ixed energy/angle mesh before proceeding to the next sublayer. 

The efficacy of this approach, i.e. , the excellent agreement with Monte 
Carlo calculations, is seen in Fig. 3. The BETA-II calculation used sepa- 
rate runs for the high {>10 MeV) aud low (<10 MeV) energy portions of the 
spectrum. The SHIELD run required about one minute of Uni vac 1108 time, 
while the PETA-II calculation required over 30 min. 

Other SHIELD capabilities include: 

(1) Energy-dependent flux output. 

(2) Angular flux output. 

(3) Multiple response functions. 

(4) Punched card output kernels for SIGMA. 

(5) Proton/heavy charged particle transport. 

^ (6) Secondary bremsstrahlung kernels. 

(7) Angular incidence, including monodirectional. 
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C. MINIMUM PATH AND ANGULAR TRANSMISSION KERNELS 



The major deficiency of slant path kernels is that they are correct 


j 

1 


only it material distributions around a receptor are actually spherically 
symmetric or if the charged particles are not deflected during tiansport. 

Figure -i indicates the error introduced by using slant path kernels 
with the Jupiter environment (Ref. 1) for simple uniform thickness slab 
geometry. In view of the uncon.rervative results (underestimates) and of the 
many spacecraft volumes which have a slablike geometry (e. g. , points just 
inside vehicle skins), an alternate one-dimensional kernel, ^ (r , flU, 

was generated by SHIELD and implemented in SIGMA. In tliis kernel, 
x^i^ (r , fl) is the minimum mass thickness path at the dose point T in the 
direction e. g. , for slab geometry, fl) is the slab thickness, 

regardless of the relative direction between O and the slab normal. This 
kernel is assumed to have an a-igular dependence of the form 

= P.pher.'H lco» 

where z = x . (r, fl), D , (z) is the dose at the center of a spherical 

shell of thickness 7 , 4> is the incident angle, and c(z) is a depth-dependent 
exponent. Bv requiri ig that this kernel correctly predict both center-of- 
sphert and slao geometry results, c(z) is simply 




and is obtained directly from SHIELD calculations (see Appendix B). As 
seen in Fig. 4, this kernel Is exact for slab geometries. 

Using the SANDYL computer program, TRW reported (Ref. 7) addi- 
tional verification of SHIELD-generated kernt .sf. Typical differential fluence 
comparisons are shown in Fig. 5 for transmission through a 2 cm slab of 
aluminum. 
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IV, COMPLEX GEOMETRY MODELING 

A description of the distribution of spacecraft materials (nriodel) is 
required by SIGMA for radiation level calculations. The material distribu- 
tion around specific dose points is obtained by ray tracing outward from 
these points.^ The ray tr rcing methods of the FASTER III (Ref. 8)/ 

BETA U/SIGMA programs required explicit descriptions of all volumes in a 
geometry, even if void. This requirement was removed because of lengthy 
input required to describe the man’ voids of complex shape in the space- 
craft. The spacecraft models used in the kernel analysis are fully compati- 
ble with Monte Carlo anal’*sis methods. 

The SIGMA and BETA-II program files both utilize a series of user- 
oriented input /outpvt data processors. The processors described below 
were used extensively. 

A. SURFACE AND REGION PROCES 'RS 

Some parts of the spacecraft wer scribed by the surface /region 
methods of the FASTER-III and BETA- rograms. Each region is defined 
by specifying the surfaces which form ii. boundaries, where each surface 
is defined by a general quadratic equation: 

. Z I Z 

G(x, y, z) - + ajX + a^y + a^z + a^x + a^y + a^z 

+ a.^xy -t agyz + a^zx =0 

where the a, are the surface coefficients, 
k 

Material regions are then defined by specifying the particular surfaces 
(simple planes, cones, cylinders, or spheres) ..Tat bound that region. 
Regions of the MJS'77 spacecraft described by this method include the pro- 
pellant tank and the high-gain antenna. 


Aluminum was used as the reference dose material. 
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The method of material distribution description givsn above was 
: simplified for use by the MJS'77 project. 

B. ELECTRONICS BAYS PROCESSORS 

The MJS'77 spacecraft has ten electronics bays of similar gecmetry. 
An input processor was written to generate the surfaces and regions com- 
prising each bay. Input parameters include: 

(1) Nurtiber of bays. 

(21 Thickness and materials of each bay wall. 

(3) Total w-.ipht and material of the bay interior. 

In geiieral, alt bays except the one of specific interest were described 
by a smeared density which conserved mass, volume and shape, 

C. BOARDS PROCESSOR 

Several electronics bays had interiors containing a series of parallel 
electronics boards. .An input processor was written to accept a simple 
input for parallel board geometries all having common transverse 
boundaries. 

D. DESIGN PROCESSOR 

Other electronics bays and all of the science instruments required 
more comiplicated geometric cescription». They were composed of many 
parts, each with a different, but simple, geometry. Again, a special input 
processor was developed to generate the surfaces and regions required to 
describe these parts. Recognized sl;ape8 include plates, cylinders, annuli, 
spheres, and truncated cones. This processor included error testing for 
overlap of legions. 

E. ROTATE AND TRANSLATE PROCESSORS 

The interior of each bay was described in a bay -centered coordinate 
system. This description was then rotated and translated, by the program, 
to the appropriate spacecraft coordinates. This procedure is schematically 
indicated in Fig. 6. 
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This same rotation and translation capability was used for detector 
points: i. e. , detectors were specified in the bay-centered coordinate system 
and then moved to the appropriate bay. 

F. CAMERA AND PICTURE PROCESSORS 

Two plotting routines were used to facilitate geometry checkout. One 
routine (PICTURE) generates printouts of geometry cross sections (Fig. 7). 
The second picture routine uses the CRT plot capabilities of the CAMERA 
program (Ref. 9 and Fig. 8). 

V. THREE-DIMENSIONAL CALCULATIONS 

A. SIGKfA PROGRAM 

Radiation levels are calculated in the SIGMA program by numerical 
integration over solid angle of the one-dimensional attenuation kernels. For 
example, the dose D(r) is given by 


D(r) = D(x(r,n))dn= ^y^K.D.(x(r,n)) An. . 

^4ir i= 1 

where D(x(r, fl)) is the dose that would be received at the dos.j point if the 
mrss thickness x(r,ii) of materials encountered in the differential solid 
angle dfl about the direction were spherically symmetric about r. The con- 
stant K. is associated with the numerical integration scheme (e. g. , K. = 1 
for midpoint integration). Typically, for MJS'77, the polar and azimuthal 
angles were each segmented into 26 divisions, making an angular integra- 
tion grid of 26 X 26 = 676 solid angle sectors (ray traces). Geometry 
mockup capabilities available in ray tracing programs such as SIGMA permit 
very accurate representation cf mass distributions. SIGMA accepts multiple 
radiation-type kernels (electron dose, proton fluence, etc. ). Ray tracing 
about a detector point is performed only once; as each ray trace is per- 
formed, the contribution to every radiation type is obtained. In particular, 

^11 SIGMA runs output both slant path (center -of- sphere) and minimum path 
(slab) results. 
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SIGMA outputs include mass path distribution (both slant path and 
minimum path) and sensitivity of radiation levels with respect to shielding 
added to particular (u ser- specified) surfaces. This shield sensitivity output 
includes recognition of unique shield crossing combinations, e. g. , none, 
one, or combinations of two or more, with corresponding output for the 
variation of the radiation level when tl ese shield thickness are varied. 
Typical dose sensitivity data for one dose point is shovm in Fig. 9. 

B. SHIELD OPTIMIZATION 

SIGMA obtains shield- sensitivity data for one or more detector points 
and optionally saves the data on a psrmanenv file. These data are then avail- 
able for shield optimization calculations. 

Optimization calculations can be performed for multiple dose points, 
shields, and criteria. The user specifies the criteria upper limits and how 
each criterion is formed from the individual radiation kernels calculated by 
SIGMA, For example, 

Dose criterion = (1) x (minimum path electron dose kernel) 

+ (1) X (slant path proton dose kernel) 

+ (0) X (all other kernels). 

The user also specifies the geometry and minimum and maximum 
shield thickness at each candidate shield location. Candidate shields are 
those specified for the original SIGMA calculation plus a unit shield for 
each point detector. Unit shields, which are spherical shields centered 
at the dose point, are not obtained if the maximum unit shield thickness is 
specified as zero. 

; The optimum shield configuration is calculated by iteratively incre- 
menting shield thickness until all criteria are met. On each iteration, the 
change in the jth radiation level is calculated separately for an increment 
of At: in the ith shie,ld; i. e. , 

AL.. = L.(t + At.) - L.(t) 

J J 
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is the change in the jth radiation level due to changing the thickness 
t = t (tj, t^f • • • tjj) in the ith shield only. 

The corresponding change in the total shield weight is calculated as 


AW. = W.(t + At.) - W.(t) 
1 1 1 1 


using weight equations for slab, cylindrical, and/or spherical geometries, 
either isolated or nested. 


A combined relative shield worth is calculated as 




Q: = 


AW. 


where the summation is over different radiation level criteria and L. is the 

J 

jth criterion. That shield for which is most negative is changed by a 
thickness At. and the process is repeated until all criteria are met. 

The results obtained for a single detector point and a single criterion 
are shown in Fig. 10. Because the optimization uses interpolations of tabu- 
lated sensitivity data, a SIGMA calculation for the optimized shield configura- 
tion agrees with the optimization output to within a few percent. 


VI. 


CONCLUSIONS 


Ray tracing (sectoring) transport programs like SIGMA do have their 
limitations. The errors introduced by using one-dimensional attenuation 
kernels and by assuming that electrons do not scatter from one solid angle 
sector to another may be significant. Unfortunately, no experimental 
measurements that can be used for direct quantitative assessment of the 
accuracy of these programs have been made for three-dimensional shield 
configurations. Nevertheless, compared to Monte Carlo programs, the 
computer programs described In this report are fast, convenient, versatile, 
and inexpensive. Together, they represent a necessary capability for any 
project where space radiation shielding engineering is an essential discipline. 
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Fig. 1, h/lariner Jupiter/Saturn 1977 spacecraft. 



Fig. 2. Electron and proton radiation transport software flow 
diagram. 
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Fig. 5. SHIELiD/SANDYL. comparison for 
tranf.Tjission through a Z-cm slab of 
alumiuium. 
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Fig. 6, Schematic iUustration of the 
: sequence of rotations and translations 
used ir. SIGMA to position an elec- 
tronics bay at its proper location in 
the spacecraft. 
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Fig. 8. Compufer plot of the photopolarimeter subsystem. The geometry model 
consists of 108 quadratic surfaces bounding 56 material regions. 
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Fig. 9. Shield sensitivities for the Fig. 10. Shield weight opiimiza- 

Canopus tracker for one dose point. tion for the Canopus tracker for 

one dose point. 
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APPENDIX A ‘ 

SHIELD, NUMERICAL METHOD 

The SHIELD code calculates charged particle transport in 
one-dimensional geometries. T^~e following definitions are used in the 
numerical integration. . 

L DIRECTION VECTORS 


Let n denote the original particle direction 


= (i cos 0 + j sin 0) •yl - p. + kp 


0 = azimuth 


p = polar angle cosine 


(A-1) 


and Ti' the deflected direction 


f0* = azimuth 


Cl' = (i cos 0' + j sin 0') vl - p'^ + kp' (Ar2) 

p' = polar angle cosins 


Then the cosine of the deflection c.nalysis 


p = n • O' = J\ - p^ Jl - p'^ cos w+ pp’ , w = 0' - 0 (A- 3) 


II. LEGENDRE POLYNOMIALS 


By definition 


= 1 . K 


(p) = i [(2£ - 1) Pj.j(p) - (e - I) P^.z^mI] (A-5i 


P.(n) = 


1 + ^^l-] 
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j: 




(A-7) 


The addition forr.;ala for azimuthal rymmetry is equivalent to 


p^(fi-n') = p^(^^) - Pj(k) p^(k') 


(A-8) 


For forward directions, the following is defined: 


1 - P^(fx) = 


(A-9) 


It follows that 


Bq(|i) = 0 . Bj(V) = » 


(A-10) 


and from (A-9) and (A-5), 


= }1(2« » 1) [l > • (f - 1) 


It can be shown by induction, using (A-ll), that 


B,,. = .. = tllfU 


(A-12) 


Therefore, 


lim j^l - (t - (x) 


(A-13) 


The expansion of a general function f(p) uses the coefficients f^, where 




f = 2tr / f(;- p^{p) dp 


(A-14) 
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1. e. 


£ = 0 


(A-I5) 


Finally, if for any function g (fi) 


g(fi) = g - --- (a delta function), then g; = g for all I (A- 16) 




III. SCATTERING CROSS SECTION’ 


The charged particle deflection cross section iia<> the form 


d<r 


(I - M + n) 


(A-17) 


where n is the screening angle, and ti « 1 => strong peaking at p = 1, 
It follows that 


/■‘si 

L 


«r = 2n / (p) dp 


(A- I 8) 


is the total, cross section and the Legendre expansion coefficients are 


<Tg = 2w 




(p) dp 


(A-19) 


in particular - (t. Thus, the scattering cross section can be repre- 
sented by 


/ 


d<T _ . 

dii - 


(A-ZO) 


;=o 
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By definition (for later use), 




f % < 1*1 [' - '“>] 


dji 


(A-21) 


or equivalently. 


j 


(^) B^(iJi) (1 - k) d|x 


(A-22) 


With strong forward peaking, most of the integral comes from p « 1, so that 




or* =s 2ir B^(p = I ^ 


(A-23) 


or 


# ill + 1) , 
^ 2 ^ 


(A- 24) 


where 


■ ”/> 


(p) (1 - p) dp 


(A-25) 


Equality in Eq. (A-24) is used for er* rather than the exact value from 
Eq, (A-21) or (A-22). 
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IV. TRANSPORT EQUATION 

The transport equation is wri*^ten as a function of particle distance 
trz.versed, where 




3 


4>(s, p) = the angular flux at distance s 
The Legendre moments of the flux are 


(A-26) 


iS-i 

& 


4^ (s) = 2tr 

m 




1a) 


(A-27) 


i 

u 


Therefore, the angular flux can be represented by 


<t'(s 


• ■‘I ' V,'»> Pm''*' 


m=0 


(A-2&) 


The transport equation, allowing only deflection reactions, ic 


||(s,p)= I <l>(s, p') |^(s,n' . ni d«' - (t{s) 4>(s. p) (A-29) 


/ 

J4ir 


Substituting for the scattering cross section yields 


rS 

$ 


i" . 
5T • 

^ ' 
i ■■ 

$■ 

n ■ 

I- 

it 

3 5 

t* 


84, , 


I 4(s, p') 

00 

P^(P) P£(P') 

/4w 

J = 0 


dp’d0' 
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- <r(s) 4>(s, p) 

(A- 30) 

ft 

>5 • 

CO 



^ 4tt (s) P^ (p)<{)^(s) - <r(s) 4 >(s, p) 

(A-31) 

.1 

i . 



- 

% • 



vh - 
' 


•f 

.4 

m 4 



Multiplying by P^(p), integrating ever 4w, and using 


-p. 


i: 


Im 


P^(^) P^(„) 


(A-32) 


yields 


‘ 4 

‘ t 


9<t> 


m 


OS mm m 


-[o^(s) - cTm(s)] <j>^(s) (A-33) 




or, using the definition of 




m 


ds 


(s) = 


-cr';' (s) 4> (s) ^ fm«|) (s) 

m m os m 


-ff^n^s) (A-34) 


A simple integral yields the solution 


<t>,„(s) = 4>_(0) exp 
m m ■ 


/■ 


<r„(s’) ds' 
m 


(A-35) 


where the <}> (0) are the Legendre moments at s = 0. Thus 

m . 




<|)j^(0) = 2ir / <j)(0, p) P^(p) dp 


(A-36) 


Equations (A-35) and (A-36), and Eq. (A-21) for o*, constitute the 
Gaudsmit-Saunder sor method. The series is truncated at order L and a 
delta function component is assumed. Therefore, 


<l>(8, p) = ${S, p) + «I>l(8) 


(A-37) 
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The moments of this equation are 

<l>_(s) = 5_(s) + 6. (s) (A-38) 

m m j_> 

Therefore, 

= 4>^(s) - <J>T (s) (A-39) 

m m 

and 

L-1 

<t.(s, ji) = - <j>j^(s)] P^(^) + ,j:^(s) Miy-U (A-40) 

m= 0 

V. TRANSPORT THROUGH MATERIALS 


The above definitions are used for the charged particle transport in 

the following numerical scheme. First, the shield material thickness is 

divided into many small layers; each layer thickness is much less than the 

range of the lowest energy electrons. Second, the range of polar angle 

cosine relative to the layer normal is divided into many intervals; i. e. , 

interval i is u. < u < u. , , where the u.. 's are the interval boundaries. The 
n ^ '^i+l 1 

transport problem is then solved by repetitively calculating the transmitted 
energy-angle distribution for a layer and using that distribution as the inci- 
dent source for the next layer. 


In particular, let denote the number of incident p?,rticles in energy 
grouo j and angle interval i. Let Ej. denote the average energy of the 
particles. The problem is a calculation of W^. and e!. or the other side of a 
layer of thickness t; i. e. , for each energy group and a*, ^le interval. 


^i ” T ^^i ^ ^i+P ’ average polar angle cosine 
s. = — , the path length across the layer 
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E.. 




E.. - I ^ (s’) 

Jl f ds ' ' 


Jo 


ds' , the transmitted energy 


n 

'ji = J ”j 
•'0 


(s') ds' , b(s) from (A-25) 


<|>ji = exp 


|-L(L+ « the fraction of particles 

L J mitted into the same an 


:s trans- 
mitted into the same angle 
interval. 


The energy-angle distribution of deflected particles is 


W:’.(E', p', 0') 


where 


W.. 

Jl 


L-I 


P (p) 

m 


rn-.= 0 


2m + 1 . AW y . . 

'^jim - tij) m 


6(E- - Ej.) 


A.. = exp I - m (m 

Jim I 


|j- m (m + 1 ) I" L(L + 1 ) • 


Integrating over the initial and final directions and adding on the undeflecte(<. 
component yields the energy distribution of transmission into angle interval i': 


fL-1 


w!',..(E') = W..[y^ A.. + <>?. 6.,.) 6(E- - E’.'.) , 

1 Jl ji|Z^ 2 jii" fii-i-i ■ >^1 J* * ‘I 


1 m=0 


where C. is the integral of the Legendre polynomial over the angle 


im 


interval. Integration over final energy yields 


i j ./group j' 
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: APPENDIX B 
MINIMUM PATH KERNEL 

For isotropic incidence on both sides of a slab. 


D(z) 


’ jt 


Zn cos 0 = 1 


cos 0=0 


Td . (z)l (co3 0)^^^^”^ d(cos 0) d0 
I sphere I 


(B-1) 


V ^ 


Td 

I sphere 

Jo 


(z) dp , p = cos 0 


^sphere^^^ 




^sphere^^^ 

c(z) 


Ps£her^ 
/Psphere<^ \ 
V ^slab(^> / 


D(z) = 


°slab<^) 


For isotropic incidence on a sphere, p = 1 and 


rZrt ^ cos 0=1 

f [“sphered] 

/O Jcos 0=0 


D(z) 


(z)l d(cos 0) dO 


so that, from (B- 1) 


D(z) = D , (z) 

sphere 
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